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Abstract. Dust formation in brown dwarf atmospheres is studied by utilising a model for driven turbulence in the 
mesoscopic scale regime. We apply a pseudo-spectral method where waves are created and superimposed within a 
limited wavenumber interval. The turbulent kinetic energy distribution follows the Kolmogoroff spectrum which 
is assumed to be the most likely value. Such superimposed, stochastic waves may occur in a convectively active 
environment. They cause nucleation fronts and nucleation events and thereby initiate the dust formation process 
which continues until all condensible material is consumed. Small disturbances are found to have a large impact on 
the dust forming system. An initially dust-hostile region, which may originally be optically thin, becomes optically 
thick in a patchy way showing considerable variations in the dust properties during the formation process. The 
dust appears in lanes and curls as a result of the interaction with waves, i. e. turbulence, which form larger and 
larger structures with time. Aiming on a physical understanding of the variability of brown dwarfs, related to 
structure formation in substellar atmospheres, we work out first necessary criteria for small-scale closure models 
to be applied in macroscopic simulations of dust forming astrophysical systems. 



1. Introduction 

Substellar objects like brown dwarfs and (extrasolar) plan- 
ets are largely - but not entirely - convective with con- 
siderable overshoot in the upper atmosphere. They also 
provide excellent conditions for the gas phase transition 
to solids or liquids (henceforth called dust). The under- 
standing of such object's atmospheres requires therefore 
the modelling of convection and, because the inertia of 
the fluid is larger than its friction (Re::^ 10^), turbulent 
dust formation must be considered. Hence, substellar at- 
mospheres involve various scale regimes with each being 
dominated by possibly different physical (e.g. streams, 
waves, precipitation) and chemical processes (e.g. com- 
bustion, dust formation, coagulation). 

The large-scale structure of compact, substellar atmo- 
spheres is characterised by local and global convective mo- 
tions (e. g. thunderstorms and monsoon like winds) and 
- simultaneously - by the gravitational settling of the 
dust. This scale regime has widely been investigated by 
ID static, frequency dependent atmosphere calculations 
applying mixing length theory. The presence of dust in 
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form of several homogeneous constituents has been mod- 
elled by applying local stability and time scale arguments 
(Rossow 1978, Borrows et al. 1997, Saeger & Sasselow 
2000, Ackermann & Marley 2001, AUard et al. 2001, Tsuji 
2002, Cooper et al. 2003). Recently, Woitkc & HcUing 
(2003a; Paper II) have proposed a consistent treatment 
of nucleation, growth, evaporation and gravitational set- 
tling of heterogeneous dust particles, which has been ap- 
plied for the first time to stellar atmosphere models in 
(Woitke & HeUing 2004; Paper III). 

Modelling macroscopic scales, Ludwig et al. (2002) 
have presented the first 3D simulations of M-dwarf con- 
vective atmospheres, applying the Large Eddy Simulation 
(LES) technique based on the experience of Nordlund 
& Stein with the solar convection. While the largest 
scales, which contain most of the energy and energise a 
small-scale turbulent fluid field, are computationally re- 
solved, smaller scales are modelled by a hyperdiffusion 
which prevents the energy to accumulate in the smallest 
scales, thereby smoothing out all small-scale structures. 
Numerically, it stabilises the flow by filtering out sound 
and fast mode waves (see e. g. Gaunt & Korpi 2001). 

Coming from the opposite site of the turbulent energy 
cascade, Helling et al. (2001; Paper I) have investigated 
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the dust formation process in the small, microscopic scale 
regime [Iref <C Hp) by direct simulations of acoustic wave 
interactions. These investigations of the dense, initially 
dust-hostile layers in brown dwarf atmospheres have re- 
vealed a feedback loop which characterises the dust forma- 
tion process: 

The interaction of small-scale perturbations of the 
fluid field can cause a short-term temperature de- 
crease low enough to initiate dust nucleation. The 
seed particles grow until they reach a size where the 
dust opacity is large enough to re-enforce radiative 
cooling which causes the temperature to decrease 
again below the nucleation threshold. Dust nucle- 
ation is henceforth re- initiated which results in a fur- 
ther intensified radiative cooling. The nucleation rate 
and consequently also the amount of dust particles 
increase further. This run-away process is stopped 
if either the radiative equilibrium temperature of the 
gas is reached or all condcnsiblc material is con- 
sumed. Meanwhile, the seed particles have grown to 
macroscopic sizes (^m). 

Based on this knowledge, we extent our studies to larger 
and larger spatial scales aiming finally at the simulation of 
dust formation in the macroscopic scale regime, i.e., the 
complete atmosphere. The next step is therefore to study 
the mesoscopic scale regime (Zief < Hp) where we model 
driven turbulence by stochastically superimposed waves 
in the incrtial Kolmogoroff range and study the response 
of the dust complex. Necessary criteria are derived for 
a small-scale closure model to be applied in large scale 
simulations of dust forming systems. 

The aim of such a scale- wise investigation is to under- 
stand the major physical mechanisms which are respon- 
sible for the structure formation in the atmospheres of 
substellar objects and to provide the necessary informa- 
tions for building an appropriate sub-grid model needed 
to solve the closure problem inherent to any macroscopic 
turbulence simulation (not only) of dust forming media 
(see also Canute 1997a, 2000). The challenge is that only 
the largest scales of the turbulence cascade are compa- 
rable with real astrophysical observations but structure 
formation is usually seeded on the smallest scales espe- 
cially if it is correlated with chemical processes. From the 
theoretical point of view, turbulence in thin atmospheric 
layers may be of quasi-2D-nature (Cho & Polvani 1996, 
Menou et al. 2003) . The 2D turbulence is characterised by 
an inverse energy cascade (transfer from small to large 
scales), contrary to 3D turbulence, which makes a scale- 
wise investigation even more urgent. 

Various model approaches have been carried out to 
simulate and to study turbulence in different astrophys- 
ical scale regimes. For example, thermonuclear flames in 
type la supernovae have been studied on small scales by 
Ropke, Niemeyer & Hillebrandt (2003), investigating the 
Landau-Darrieus instability in 2D simulations which is re- 
sponsible for the formation of a cellular structure of the 
burning front. Reinecke, Hillebrandt & Niemeyer (2002) 



have performed large-scale calculations in order to model 
supernovae explosions on scales of the stellar radius. Mac 
Low (1999; see also Mac Low & Klessen 2004), for in- 
stance, have set up isothermal initial velocity perturba- 
tions with an initial power spectrum of developed turbu- 
lence in the Fourier space and initially constant density. 
They model decaying turbulence on the small scales of 
the inertial subrange. Smith, Mac Low & Heitsch (2000) 
use a similar approach but study the effect of driven tur- 
bulence in the same scale regime of star-forming clouds. 
A stationary but stochastic velocity field was applied by 
Wallin, Watson & Wylf (1998) in order to perform radia- 
tive transfer calculations of Maser spectra of a sub-parsec 
disk of a massive black whole. A fundamental theoreti- 
cal investigation of the methods of driven turbulence is 
provided in (Eswaran & Pope 1987, 1988). A different ap- 
proach of turbulence and convective modelling is followed 
by Canuto (e. g. 1997b) who treats turbulent convection 
by a Reynold separation ansatz where decomposed quanti- 
ties (background field -|- fluctuations) are introduced into 
the model equations. 

In this paper, we present a model for driven turbulence 
which allows us to study the onset of dust formation un- 
der strongly fluctuating hydro- and thermodynamic con- 
ditions in the mesoscopic scale regime. We thereby intend 
to model a constantly occurring energy input from the 
convectively active zone. Section Estates the model prob- 
lem and the characteristic numbers of our astrophysical 
problem. The turbulence model is outlined in Sect. 12.21 In 
Sect. 121 the results are presented for ID simulations and an 
illustrative 2D example. Section^jcontains the discussion, 
and Sectional the conclusions. 

2. The model 

The model equations for a compact substellar atmosphere 
arc summarised. The model is threefold: i) hydro- and 
thermodynamics {Complex A), ii) chemistry and dust for- 
mation (Complex B] both Sect. I^TTjl . and iii) turbulence 
(Sect. I2.2|l . Our model philosophy is to study the dust 
forming system in different scale regimes in order to iden- 
tify major mechanisms which might be responsible for 
cloud formation and possible variability in substellar at- 
mospheres. The approach is based on dimensionless equa- 
tions such that their solution is characterised by a set of 
characteristic numbers. 

2.1. The model for a compressible, dust- forming gas 

The complete set of model equations was outlined in detail 
in Paper I. Only a short summary is therefore give here. 

Complex A: The hydro- and thermodynamics are de- 
scribed following the classical approach for an inviscid, 
compressible fluid; (Eqs.(l)-(5)) in Paper I^. 

^ Note that Eq. (2) in Paper I should be corrected as 
ipv)t + V • (p« o «) = --iyVP - ^pg. 
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Complex B: The chemistry and dust formation. The 
dust formation is a two step process - nucleation and 
growth (Gail et al. 1984; Gail & Scdlmayr 1988) - and 
depends through the amount of condcnsible species on 
the local density and chemical composition of the gas 
which are determined by Complex A. The nucleation rate 
J*, which is strongly temperature-dependent, is calcu- 
lated from Eq. (17) in Paper I applying the modified 
classical nucleation theory^ (Gail et al. 1984). The dust 
growth is described by the momentum method developed 
by (Gail & Sedlmayr 1986, 1988; Dominik et al. 1993) in 
combination with the differential equations describing the 
element conservation (Eqs.(6)-(8) in Paper I). 

Our model of dust formation considers a prototype- 
like phase transition (gas solid) which is trig- 
gered by the nucleation of homogeneous (Ti02)jv- 
clusters (see Jeong et al. 1998, Gail & Scdlmayr 1998, 
Jeong et al. 2003). The formation of the dust particles is 
completed by the growth of a heterogeneous mantle which 
is assumed to be arbitrarily stable. The most abundant el- 
ements after H, C, O, and N in a solar composition gas are 
Mg and Si followed by Fe, S, Al, . . ., Ti, . . . Zr. Therefore, 
the main component of the dust mantle can be expected 
to be some kind of silicate with a Mg/Si/0 mixture plus 
some impurities. Since the focus of our work concerns the 
initiation of the dust formation in hostile turbulent envi- 
ronments rather than a detailed description of the growth 
process, evaporation and drift are neglected. Therewith, 
the maximum effects regarding the amount of dust formed 
in brown dwarf atmospheres is studied. 

The most abundant Si bearing species in the gas 
phase under conditions of chemical equilibrium is the SiO 
molecule. Therefore, the collision rate with SiO can be 
expected to limit the growth of various silicate mate- 
rials (like SiOa, Mg2:rFe2(i_:,)Si04 and Mg^Fei-^^SiOa) 
rather than the collision rate with the nominal molecules 
(monomers) which are usually much less stable and hence 
barely present in the gas phase (Gail & Sedlmayr 1999). 
Consequently, SiO is identified as the key species for the 



^ Classical nucleation theory has often been criticised (e. g. 
Michael et al. 2003) but no other consistent and in hydrody- 
namic simulations applicable theory for phase transitions was 
proposed so far. The most accurate way is the solution of the 
complete chemical rate net work for which, however, the neces- 
sary data are simply not available. The conceptional weakness 
of the classical nucleation theory is mainly the use of the bulk 
surface tension a to express the binding defects on the surface 
of small clusters. However, Gail et al. (1984) have proposed the 
modified classical nucleation theory where the bulk surface ten- 
sion is not used. Instead, thermodynamic data for individual 
clusters are adopted which are provided by extensive quan- 
tum mechanical calculations (see e. g. Jeong et al. 1998, 2000; 
Chang et al. 1998, 2001; Patzer et al. 2002; John 2003). In con- 
trast, because the calculation of high-quality thermodynamic 
data is a challenging problem, most people base their dust for- 
mation considerations simply on stability arguments (see e. g. 
Tielens et al. 1998) which is, however, a necessary but not a 
sufficient condition for phase transitions. 



description of the growth of the dust mantles in our model 
(compare Paper II). In order to prevent an overproduction 
of seed particles, we furthermore include Ti02 also as ad- 
ditional growth species, which leads to a quick consump- 
tion of Ti from the gas phase as soon as relevant amounts 
of dust are present. 

2.1.1. Characteristic numbers and scale analysis 

The use of dimensionless equations (Eqs. (l)-(8) in 
Paper I) provides the possibihty to characterise the sys- 
tems behaviour by non-dimensional numbers. The related 
estimations of typical time and length scales are sum- 
marised in Table lXdl f ApDcndix [X|l for which the reference 
values only need to be known by orders of magnitude. 
It would, however, be very difficult to adopt an unique 
representation of the reference values from recent brown 
dwarfs model atmosphere calculation because of the differ- 
ences among the different groups (compare Fig. I10|l . The 
agreement is nevertheless good enough that we consider 
them as typical, classical hydrostatic brown dwarf model 
atmosphere which lead our choice of reference values in 
Table Um 

Complex A: - Assuming the typical turbulence velocity to 
be of the order of one tenth of the velocity of sound leads to 
a Mach number Al « 0(0.1). This choice has been guided 
by the results of Ludwig et al. (2002) who derived a max- 
imum vertical velocity of ©(IC^cms"^) ~ c^/lO cms~^ 
which is about the same order of magnitude like the con- 
vective velocities derived from the mixing length theory 
(MLT)^. This value of a large-scale velocity does presently 
determine the energy dissipation rate (Eq. |3J) of the tur- 
bulence model applied fSect. in this paper. 

- The Froude number is Fr ~ 0(10^^ . . . 10^^) for 
a mesoscopic reference lengths l^d < Hp. Therefore, the 
pressure gradient and the gravity are now of almost com- 
parable importance. However, gravity will gain consid- 
erable influence on the hydrodynamics only for scales 
regimes ^ Hp {Ij-d = Hp ^ Fr ~ M). The analy- 
sis of the characteristic combined drift number performed 
in Paper III (see Tables 2, 4 therein) has shown that the 
drift term in the dust moment equations is merely influ- 
enced by the gravity and the bulk density of the grains. 
We therefore assume also in the mesoscopic scale regime 

Model atmosphere calculations for Brown Dwarfs using 
MLT provide a typical (static) convective velocity uj^^nj ~ 
0(10^ . . . lO^cm/s) Cs/1000 . . . Cs/lOOcm/s which strongly 
contradicts the value used to model an additional spectral line 
broadening component, the so-called micro-turbulence veloc- 
ity Umicro ~ Ikm/s. While ii^nJ is needed to calculate an ad- 
equate temperature structure of the inner atmosphere, Umicro 
is needed for a best fit of the observed spectra. This dilemma 
can, however, not be solved without either a consistent con- 
vective theory or a direct numerical simulation (DNS) of the 
turbulence and the convective energy transfer which influences 
the local and the global temperature and velocity fields and 
thereby the observed spectral lines and their broadening. 
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position coupling between dust and gas which seems rea- 
sonable due to the almost equal importance of the source 
terms in the equation of motion. 

- The estimate of the Reynolds number, Re = 
10^ . . . 10^, for a brown dwarf atmosphere situation in the 
mesoscopic scale regime indicates that the viscosity of the 
gas is too small to damp hydrodynamical perturbations 
on the largest scale to be considered, Zmcso, and a tur- 
bulent hydrodynamic field can be expected. Re has in- 
creased by about one order of magnitude compared to the 
microscopic scale regime (compare Table 1 in Paper I). 
Therefore, the viscosity of the gas decreases for mesoscopic 
scale effects in comparison to the microscopic regime. This 
is correct since for l^-d ^ rj ^ Re = 1 (77 - Kolmogoroff 
dissipation scale) and viscosity dissipates all the turbulent 
kinetic energy of the fluid. 

- In radiatively influenced environments the charac- 
teristic number for the radiative heating / cooling, Rd = 
4K;ref cTf^f ' ■ The systems scaling influences Rd by the 
reference time trei which can increases with increasing spa- 
tial scales. 



Complex B: The scaling of the dust moment equations 
provides two Damkohler numbers for dust nuclcation, 
Da^^^, and dust growth. Da^ , and characteristic num- 
bers for the grain size distribution, the Scdlmayr number 
Scj (j - order of dust moments, j e N). Element conserva- 
tion is characterised by El, the clement consumption num- 
ber. Scj and El are not influenced by the scaling of the 
system (compare Table TKA^ but the two dust Damkohler 
numbers, Da^^'^ and Da^^ , increase with increasing time 
scale. 

The analysis of the characteristic numbers shows that 
the governing equations of our model problem arc still 
those of an inviscid, compressible fluid which are coupled 
to stiff dust moment equations and an almost singular 
radiative energy relaxation if dust is present. The dust 
equations become even more stiff than in the microscopic 
regime which caused severe numerical difficulties in solv- 
ing the energy equation which is coupled to the dust com- 
plex by the absorption coefficient k (see Sect l2.3.T|) . The 
dominant interactions occur in the energy equation and in 
the dust moment equations (Eqs. (4), (6), (7) in Paper I) 
also in the mesoscopic scale regime. 



2.2. The model for compressible, driven turbulence 

A turbulent fluid field is determined by the stochas- 
tic character of the hydrodynamic and thermodynamic 
quantities due to possible interaction of different scales 
which are represented by inverse wavenumbers in our 
model. We have constructed a pseudo-spectral method 
where randomly interacting waves are generated inside a 
wavenumber interval [/cmin, ^max] on an equidistant grid of 
N wavenumbers fcj in the Fourier space. The wavenumber 
interval is part of the inertia! subrange of the turbulent 
energy cascade (Eq. EJ . 



A disturbance 5a{x,f) is added to a homogeneous 
background field ao{x,t) such that for a suitable variable 



y.{x, t) = ao(a;, t) + 5a{x, t). 



(1) 



with X the spatial vector and t the time coordinate. The 
present model for driven, compressible turbulence com- 
prises a stochastic, dust-free velocity, pressure and entropy 
field, i.e. a{x,t)e{u{x,t),p{x,t),S{x,t)}. 

Stochastic distribution of velocity amplitudes 5u{x,t): An 
arbitrary scale - represented by a wavenumber interval 
k . . .k + dk [k = \\k\\) - inside the inertial range of devel- 
oped turbulence contains the energy per mass e{k) dk in 
3D, where 

e{k)dk = CKe^/^k~^/^. (2) 

£ is the energy dissipation rate [cm^/s'^] and Ck ~ 1.5 
is the (dimensionless) Kolmogoroff constant (see Dubois, 
Janbcrteau, Teman 1999, p. 51). Kolmogoroff derived 
this first order description of the energy spectrum for 
the inertial range assuming self-similarity of the corre- 
sponding scales. The energy spectrum (Eq. ^ has been 
well verified by experiments and simulations (see e. g. 
Dubois, Janberteau, Teman 1999). In the framework of 
Kolmogoroffs theory, e is constant for all scales k and time 
t {homogeneous, isotropic turbulence). 

The energy dissipation rate can be estimated if already 
one typical scale and its corresponding reference velocity 
is known because e is assumed to be constant for all scales 
(see also Sect. 12.111 Comvlex A). From dimensional argu- 
ments. 



(3) 



for instance for the largest scales of interest inside the 
inertial range, i.e. for the smallest wavenumber /cmin- 
According to Jimenez et al. (1993), Ci = 0.7. 

A wavenumber interval [/cijfcj+i] contains, according 
to Eq. Q, the turbulent kinetic energy density per mass 



E, 



turb' 



-^turb 



e{k)dk ^ -CKe'^-^[k. 



2/3rL-2/3 



k; 



-2/3i 



i+1 



(4) 



The square of the velocity amplitude, Au^(ki), is correlated 
with the turbulent kinetic energy in Fourier space by 



Au(fci) = J2Z3EI 



^turb' 



with 



h = 



ki -\- ki 



+1 



(5) 



(6) 



the mean value of k in the wavenumber interval consid- 
ered, fci are N equidistantly distributed wavenumbers in 
the Fourier space, (i= 1,...N) with N the number of 
modes. The fcj are chosen between /cmin and /cmax when 
the calculation is started and are kept constant further 
on. 
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Here, the so-called ultraviolet truncation Au{ki) = for 
ki > fcmax is applied in order to avoid the infinite energy 
problem of the classical field theories in Eq. (stated 
in Bohr 1998, p. 23). The minimum wavenumber is deter- 
mined by the largest scale l^^f, i- c. the size of test volume. 
Only wavcnumbcrs inside a sphere of radius fcmax exclud- 
ing the origin are forced (see also Ovcrholt & Pope 1998, 
p.l3). 

Assuming the ergodic hypothesis (see e. g. Frisch 
1995), the turbulent kinetic energy El^^.^^ (Eq. ^ was as- 
sumed to be the most likely value (compare e.g. Mac 
Low et al. 1999) with a stochastic fluctuation generated 
by a zero-centred Gaussian distributed random number 
Z3 > according to the Box-Miiller formula 



Z3 



\/ -1 log^i sin(7rz2) 



(7) 



z\ and Z2 are equally distributed random numbers e [0, 1). 

A cosine Fourier transformation provides the real val- 
ues of the velocity amplitude Su{x,t) in ordinary space. 



5u{x, t) = Si Au{ki) cos(feia; — cjjt + (pi)k, 



(8) 



with ki = kik and u{x,t) = u{x,t)u. Also the directions 
of fci, i.e. the direction of 6u{x,t), are chosen randomly 
according to 

fci.x = sin a cos/3 cos a = 1 — 2z4 (9) 

fci^y = sina sin/3 sin a = ^1.0 — (cosa)^ (10) 

fci,z=cosQ! /3=27rz5, (11) 



with Z4 and Z5 equally distributed random numbers. The 
ID and 2D case of Eq. ijSJ) is obtained by projection. A 
longitudinal wave results in ID. 

ipi = 27r zg is the equally distributed random phase 
shift which is chosen separately for each wavenumber. Wj 
is the angular velocity for which a dispersion relation is de- 
rived from dimensional arguments. It follows from Eq. (|3J) 
that for each scale h ~ 27r/fci the corresponding eddy 
turnover time t\ results to be 



Mi 



ti 



-1/3 



(12) 



Since per definition loi = 271 /ti the dispersion relation in 
the inertial subrange is 



c^i = (27rfcfe)i/3. 



(13) 



Stochastic distribution of pressure amplitudes Sp{x,t): 
The pressure amplitude is determined depending on the 
wavenumber of the velocity amplitude Au(fci) such that 
the compressible (sound waves) and the incompressible 
pressure limits are matched for the smallest fcmin and the 
largest fc,nax wavenumber, respectively, 



Ap(fci) 



[kn 



ki] pAu(fci)2 + [ki - fcmin] yOCsAu(fci) 



[kn 



The maximum wavenumber fcmax = 27r/(3Aa;) is de- 
termined by some factor (here 3; see Overholt & Pope 
1998 for discussion) of the spatial grid resolution Aa; (see 
Table 1X2)1 . 

A spectral decomposition (compare Eq.lHJ provides the 
real values of the pressure amplitude 5p{x, t) in ordinary 
space, 



5p{x, t) — — Ei Ap(fci) cos(fcia; — Uit + ipi). 



(15) 



(14) 



Stochastic distribution of the entropy S{x, t): The entropy 
S{x,t) is a purely thermodynamic quantity and a distri- 
bution can in principle be chosen independently from the 
distribution of the hydrodynamic quantities. In the adia- 
batic case, S{x,t) is conserved along particle trajectories. 
So far, S{x,t) has been kept constant. 

For a given S{x, t) and p{x, t) (see Eq.^l the gas tem- 
perature T(a;, t) is given by 

iogr(.,o^^^"'^) + ^^°gP^"^^)-^^°g^, (16) 

with R the ideal gas constant, cy specific heat capacity 
for constant gas volume. 

In this work, we have simulated the turbulence 
by prescribing boundary conditions (Sect. I2.3|l applying 
Eqs. ((HJ), ((T5|l . and lfTC|) . Stochastically created and su- 
perposed waves continuously enter the model volume and 
are advcctivelly transported inward by solving the model 
equations (Scct.[^^. A hydrodynamically and thermody- 
namically fiuctuating field is generated which influences 
the local dust formation because it sensitively depends on 
the local temperature and density. 

2.3. Numerics 

The fully time-dependent solution of the model equations 
has been obtained by applying a multi-dimensional hydro 
code (Smiljanovski et al. 1997) which has been extended in 
order to treat the complex of dust formation and elemen- 
tal conservation (Eqs. (l)-(8). Paper I). The hydro code 
has already been the subject of several tests and studies 
in computational science (see also Schneider et al. 1998, 
Schmidt & Klein 2002). 

Boundary conditions and turbulence driving: The Cartesian 
grid is divided in the cells of the test volume (inside) and 
the ghost cells which surround the test volume (outside). 
The state of each ghost cell is prescribed by our adia- 
batic model of driven turbulence fSect. l??^ for each time. 
Hence, the actual fluctuation amplitudes of the fluid fleld 
{5u{x,t), 6p{x,t), 6S{x,t) ST{x,t), 5{px,t)) is the 
result of a spectral composition of a number of Fourier 
modes which are determined by the Kolmogoroff spec- 
trum. The absolute level of this energy distribution func- 
tion is given by the velocity ascribed to the largest, i. e. 
the energy containing scale of the simulation (see Complex 
A in Sect. im)! . 
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Fig. 1. Time sequence of 4 time steps for a ID simulation with Tj-ct = 2100K, M — 0.1, iVk = 500( y entree A 
in Table IX!2l black/blue solid - 0.48s, grey/cyan solid - 0.68s, black/blue dash-dot - 1.12s, grey/cyan dashed-dot - 
1.5s). The first instant of time shows the superimposed waves which just enter the test volume. The later times show 
nucleation fronts {t = 0.68s) and nucleation events {t > 0.5s) occurring. 

T, p, u, Pgas are given dimensionless, the dust quantities have their physical units (J*/?i<h> in [s^^], Xhct in [cm s^^], 
(a) in [cm], rid in [cm"^], /ti in [-], /si in [-]) 



The hydro code solves the model equation in each cell 
(test volume + ghost cells) and the prescribed fluctuations 
in the ghost cells are transported into the test volume by 
the nature of the HD equations. The numerical boundary 
occurs between the ghost cells and the initially homoge- 
neous test volume and are determined by the solution of 
the Riemann problem. Material can flow into the test vol- 
ume and can leave the test volume. The solution of the 
model problem is considered inside the test volume. 



Initial conditions: The (dimensionless) initial conditions 
have been chosen as homogeneous, static, adiabatic, and 
dust free, i. e. po = 1, po = 1, mo = 0, Lq = (^ Lj = 0) 
in order to represent a (semi-)static, dust-hostile part of 



the substellar atmosphere. This allows us to study the in- 
fluence of our variable boundaries on the evolution of the 
dust complex without a possible intersection with the ini- 
tial conditions. 

2.3.1. Stiff coupling of dust and 
radiative heating / cooling 

Dust formation occurs on much shorter time scales than 
the hydrodynamic processes (see e. g. Sect. 2.2 in Paper I). 
Approaching regimes of larger and larger scales makes this 
problem more and more crucial (Sect. I2.1.1|l . Therefore, 
the dust moment and element conservation equations 
(Complex B) are solved applying an ODE solver in the 
framework of the operator splitting method assuming 
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Fig. 2. Time evolution in the cell centre (Tjcf = 2100K, M = 0.1, N]^ = 500 / entree A in Table rO|) . 

Left: During the first 6s (each 100th time step plotted). Right: Over 550 s« 9min (each 1000th time step plotted) 

The units are the same like in Fig. ^ 



T, p =const during ODE solution. In Paper I wc have used 
the CVODE solver (Cohen & Hindmarsh 2000; LLNL) 
which turned out to be insufficient for the mcsoscopic scale 
regime which we attack in the present paper. CVODE 
failed to solve our model equations after the dust had 
reached its steady state (compare Fig. 3 in Paper I). 
Therefore, it was not possible to simulate the equilib- 
rium situation of the dust complex in the mcsoscopic scale 
regime by using CVODE which in other situation has been 
very efficient. 

The LIMEX solver: The solution of the equilibrium situa- 
tion of the dust complex is essential for our investigation 
since it describes the static case of Complex B when no 
further dust formation takes place (i. e. where the source 
terms in Eqs. (6)-(8) in Paper I vanish). The reason may 
be that all available gaseous material has been consumed 
and the supersaturation rate 5 = 1 or the thermodynamic 
conditions do not allow the formation of dust. The first 
case involves an asymptotic approach of the gaseous num- 
ber density (or element abundance; see Eqs. 8 in Paper I) 
of S' = 1 which often is difficult to be solved by an ODE 
solver due to the choice of too large time steps. However, 
the asymptotic behaviour is influenced by the tempera- 
ture evolution of the gas/dust mixture which in our model 
is influenced by radiative heating/cooling (see r.h.s. of 
Eq. (3) in Paper I). Since the radiative heating/cooling 
(Qrad = Rdhi{T^^ — T'^) hcating/coofing rate) depends 
on the absorption coefficient n of the gas/dust mixture 
which strongly changes if dust forms. Consequently, the 
radiative heating/cooling rate is strongly coupled to the 



dust complex which in turn depends sensitively on the 
local temperature which is influenced by the radiative 
heating/cooling. It was therefore necessary to include also 
the radiative heating/cooling source term in the separate 
ODE treatment for which we adopted the LIMEX DAE 
solver. 

LIMEX (Dcuflhardfe Nowak 1987) is a solver for hn- 
early implicit systems of differential algebraic equations. 
It is an extrapolation method based on a linearly implicit 
Euler discretisation and is equipped with a sophisticated 
order and step-size control (Deuflhard 1983). In contrast 
to the widely used multi-step methods, e. g. OVODE, only 
linear systems of equations and no non-linear systems have 
to be solved internally. Various methods for linear system 
solution are incorporated, e. g. full and band mode, general 
sparse direct mode and iterative solution with precondi- 
tioning. The method has shown to be very efficient and 
robust in several fields of challenging applications in nu- 
merical (Nowak ct al. 1998, Ehrig et al. 1999) and astro- 
physical science (Straka 2002). 

3. Results 

The simulations presented in the following are charac- 
terised by the reference parameter set or the set of di- 
mensionless numbers given in Table E3 (Appendix El , 
and are carried out with an spatial resolution of = 500 
if not stated differently. After a detailed investigations of 
our ID models, the mean behaviour of the dust forming 
system is studied fSect. rO| which might, nevertheless, be 
an easier link to observations. The mean values do, fur- 
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thermore, provide a first insight regarding significant fea- 
tures of our dust forming system which a sub-grid model 
for a follow up large-scale simulation should reproduce. 
Turbulent fluctuations are discussed in terms of apparent 
standard deviations. Section will demonstrate the ex- 
istence of a stochastic and a deterministic dust formation 
regime in turbulent environments, beside a regime where 
dust formation is impossible, i. e. the problem of the dust 
formation window is discussed for substellar atmospheres. 

Section 13.51 will illustrate how stochastically superim- 
posed waves trigger the dust formation process in 2D. 
Large-scale (inside the mesoscopic regime) hydrodynamic 
motions seem to gather the dust in larger and larger struc- 
tures which is a result of multi-dimensionality. The ID 
simulations provide, however, the tool to gain detailed in- 
sight into the interactions of chemistry and physics for 
which multi-dimensional simulations are far to complex. 

3.1. Short term evolution 

An inviscid, astrophysical test fluid in ID with Ticf = 
2100 K and Af = 0.1 (/■ entree A Table is excited in 
the wavenumber interval [fcmin, ^max] by 500 modes, i.e. 
Nk = 500, and its short term evolution is demonstrated 
(Figs, n] 12 l.h.s.). The smallest eddy has than a size of 
■^min = 1 ui. The simulations assume a ID test volume in 
horizontal direction and gravity does therefore not influ- 
ence our ID results. 

Spatial evolution: Stochastically created waves move into 
the ID test volume from both sides {t = 0.48 s, Fig. ^ 
with a maximum velocity amplitude of ©(lO'^cm s~^) rep- 
resenting the turbulent velocity fluctuations. At some in- 
stant of time, the temperature disturbance due to inward 
moving superimposed waves is large enough that the nu- 
cleation threshold temperature is crossed locally (T < T5, 
Paper I below Eq. (21); t ~ 0.68s Fig. grey/cyan solid 
line). As the temperature disturbance penetrates into the 
test volume, a nucleation front forms which moves into 
the dust free gas of the test volume and leaves behind 
dust seeds which can grow to considerable sizes (compare 
the change of \ognd from t = 0.68s (grey/cyan solid) and 
t — 1.12 s (black dash-dot ) between x = and x = 0.16 
Fig. CI). 

The superimposed waves which enter the test vol- 
ume through its boundaries will also interact with each 
other after some time. An event-like nucleation results 
{t = 1.12s, X = 0.14 Fig. HJ. More dust is formed and 
meanwhile, the particles are large enough to re-initiate nu- 
cleation by efficient radiative cooling due to the strongly 
increased opacity (t = 1.5s Fig.Q] grey/cyan dash-dot). 

The result is a very inhomogeneously fluctuating distri- 
bution in size, number and degree of condensation of dust 
in the test volume when the dust formation dominated 
the dynamics of the system. The fluctuations arc stronger 
in the beginning of the simulations and homogenise with 



time (Fig. 121 l.h.s.). The long term behaviour will be dis- 
cussed in Sect. 13.21 

Time evolution: For better understanding of the time evo- 
lution of the hydrodynamic and dust quantities in a 
stochastically excited medium, the time evolution in the 
test volume's centre {cell centre) is depicted in Fig. [2 
(l.h.s) for the first 6s of the reference simulation with 
Nk = 500. 

The dust complex reaches a steady state after about 
t « 1.8 s in the centre of the test volume in the present sim- 
ulation, for which only one singular nucleation event has 
been responsible (3rd panel, Figs. |21 l.h.s.). Consequently, 
.hi = fsi = 1- In contrast, the hydro- and thermody- 
namic quantities (gas density, p, gas pressure, p, and ve- 
locity, u) continue to fluctuate considerably around their 
initially homogeneous values. The small variations of the 
mean particle size (a) and the number of dust particles 
are partly caused by the hydrodynamic motion in and out 
of the cell centre of the dust forming material and partly 
by the turbulent fluctuations themselves. The thermody- 
namic behaviour of the dust changes from adiabatic to 
isotherm as result of the strong radiative cooling by the 
dust. Therefore, the temperature drops and reaches the 
radiative equilibrium level (T = Tre)- Exactly the same 
qualitative behaviour was observed from Fig. 3 in Paper I. 

We conclude a different distribution of dust inside an 
initially dust free gas element: While in the centre of a 
gas element the dust formation process is completed (f = 
1 ) after it was initiated by waves which are emitted by 
its surrounding, disturbances from the boundary prevent 
the boundary layers of the gas elements to reach f — I. 
Consequently, a convectively ascending initially dust free 
cloud can be excited to form dust by waves running through 
it. Therefore, a cloud can be fully condensed much earlier 
than by any classical, static model predicted. 

3.2. Long term evolution 

The long term behaviour of our dust forming system sets 
in after the dust formation process is complete (/ = 1) 
and radiative equilibrium (T = Tre) is reached. Due to 
the strong cooling capability of the dust, only small devi- 
ations occur from the radiative equilibrium if compression 
waves occur which may be seen as colliding small-scale 
turbulence elements. The general change of the tempera- 
ture T Tre caused an increase of the density in the test 
volume (density level p > 1 Fig.Elr.h.s.) in order to main- 
tain pressure equilibrium (pressure level at p « 1, Fig. [21 
r.h.s.). 

The long term behaviour of p, p, and u are charac- 
terised by strong fluctuations constantly generated by our 
turbulence driving. In Fig. [21 (l.h.s), which depicts a higher 
time resolution of Fig. [21 (r.h.s.), single waves (turbulence 
elements) are still distinguishable of which only spikes out 
of a jungle of noise are left to observe in the long term 
behaviour in Fig. [21 (r.h.s.). 
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Fig. 3. Mean values (T^f = 2100K, M = 0.1, Ny, = 500 / entree A in Table rO)l . 

Left: Time evolution of the space means. Right: Time means as function of site over 10^ time steps. 



Comparable small fluctuations of the mean particles 
size, (a), and the number of dust particles, n^, occur over 
a long time. We recover here the 20% fluctuation which 
was already observable in Fig. [21 (l.h.s). Since the dust 
formation process is complete, these fluctuations must be 
of hydrodynamic origin, i. e. caused by the movement of 
the small-scale turbulence elements. 



3.3. The mean behaviour in space and time 

The mean behaviour of a turbulent, dust forming gas is 
studied. The space mean is the average over the test vol- 
ume at each time step 



N^=100 

100 200 300 400 5000 
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N, 



(17) 



and the time mean is the mean of each mesh cells over 
time, 

1 ^' 

{a)tix) ^ —Y.a,{x,U). (18) 

* i=0 
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Both represent the most plausible values of the quantity 
a{x,t) i) at a certain instant of time fEo. I17|) . and ii) 
at a certain site in the test volume (Eq. I18|) . The space 
means are calculated by leaving out the cells close to the 
boundary in order to exclude the fluctuations in the dust 
quantities due to inflowing dust-free material. Figure 13 
(l.h.s.) depicts the space means ( ■ )x as function of time 
t [s], and Fig. 121 (r.h.s.) depicts the time means ( ■ )t as 
function of x space [^rof]- 

Figures 13 shows that the space and the time mean val- 
ues differ considerably for the hydrodynamic quantities: 



Fig. 4. Same like Fig. (l.h.s.) but Nk = 100. 

Strong fluctuations of the space means (l.h.s.) occur as 
function of time while the time means (r.h.s.) exhibit com- 
paratively smooth variation. These fluctuations increase 
with increasing number of excitation modes which shows 
that the fluctuations are of hydrodynamic origin (see also 
Sect. im)) . 

Space means: The study of the long-term behaviour of the 
space means (l.h.s., Figs.Ol) discloses a considerable vari- 
ation of the hydrodynamic mean quantities. In contrast. 
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the dust quantities are almost constant in time after the 
dust formation process is completed. This result is a conse- 
quence of the assumed symmetry (ID), where every wave 
has to cross the whole test volume which is not the case 
in a multi-dimensional fluid field (compare Sect. 13. 5|) . 

The formation of dust causes the temperature to 
change considerably towards the radiation equilibrium 
level causing thus the density level to change (e. g. increase 
if T decreases) in order to recover the pressure equilibrium. 
Therefore, initially small perturbations in a dust forming 
system have a large effect on its overall hydrodynamic 
structure. 

The strong variation of the dust quantities during the 
beginning of the simulation is smeared out with increas- 
ing averaging time. Therefore, observing a spatially un- 
resolved dust forming system over a long time will not 
unhide the inhomogeneous behaviour though it will have 
a profound influence of the large-scale structure of any 
dust forming system, e.g., by the transition adiabatic — *■ 
isothermal behaviour, by backwarming in an substellar 
atmosphere, or by the enrichment and the depletion by 
gravitational settling. 



R =500, 



= 2100K, 
30 40 



f. -0.4 




Fig. 6. Same like Fig. El (r.h.s.) but M = 0.02 (/ entree 
B Table rO|) 



Time means: The variation of the hydrodynamic mean 
values is less strong in space (r.h.s., Figs.E)) than in time 
(l.h.s., Figs.OJ and resembles more common expectations 
for such average quantities than the hydrodynamic space 
means do. The density shift p(t = 0) p{t{T = Tre)) 
due to T ^ Tre is easier observable than for the space 
mean. 

The time mean of the nucleation rate, however, dis- 
closes the appearance of nucleation fronts and nucleation 
events: Waves which enter the test volume and already 
carry a temperature disturbance with T < Ts result in a 
nucleation front (e. g. a; < 0.025 ?ref r.h.s., Figs.O. Waves, 
i. e. turbulence elements, which interact inside the test vol- 
ume and only there create T < Ts for a short time result 
in nucleation events as the peak like (J*)t shows. 

Else, the dust quantities are constant in almost the 
whole test volume which is in agreement with their time 
averages. Deviations from these almost constant values 
occur only near the volume's boundaries since here fresh, 
uncondensed material enters. 

Viewing our test volume again as mass element in a 
convective environment which is constantly disturbed by 
wave propagation, we conclude that nucleation will take 
place everywhere in the mass element but likely with very 
much different efficiency. Since fresh, uncondensed ma- 
terial enters the mass element through open boundaries, 
nucleation can go on here only if the temperature is low 
enough. This does not cause the boundary region to con- 
tain the largest amount of dust since the dust can also 
leave the mass element if the fluid flow moves outward. 



3.3.1. Dependence on the number of modes 

Figure depicts the same calculation like Figs. El but 
carried out with different a number of modes, = 100. 
Comparing with the l.h.s. of Fig. (21 (iVfc = 500) shows that 
the variations in the hydrodynamic quantities are smaller 
but the dust quantities reach very alike mean values inde- 
pendent of Nk. 

Note, more energy is contained in the small wavenum- 
bers (= large spatial scales) since the Kolmogoroff spec- 
trum is applied to calculate the velocity disturbances in 
the Fourier space. Consequently, if a number of chosen 
modes, N^, is small, the smallest wavenumber will con- 
tain less energy than the smallest wavenumber for some 
larger number of modes N^, and, from iVjJ < follows 
i?(fci(N^)) < £;(fci(N^)). This results in the appearance 
of larger velocity and pressure peaks with increasing iVk 
according to Eqs. lO and 

The study of the long term behaviour of the Tiof = 
2100 K simulation reveals the occurrence of a long term 
pattern in p, p, and u (beat frequency oscillations) with 
a frequency t'boat ~ 100 s w 1.7 min. Comparing p, p, and 
u in Figures 01 exhibits 6 maxima at t « 50s, 150s, 250s, 
350s, 450s, 550s. This beat frequency v^cat seems inde- 
pendent on the number of modes Nk but does, however, 
not establish for an excitation with a very small number of 
modes (e. g. Nk ~ 5, not depicted here) and is smeared for 
a very large number of modes due to larger fluctuations 
around the mean values [N^ = 500, Fig.|2|). 
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Fig. 5. The space- means {a{t))x (solid/red) with the apparent standard deviations erg (dotted/dashed) as func- 

tion of time for T.^f = 2100 (A, l.h.s.), and T^f = 1900 (C, r.h.s.); sec Table rOl + c^S^-il^) - dotted/blue; 

{a{t))a: - o-g - dashed/green). 



3.3.2. Apparent standard deviation 

Deviations from the most plausible values, the mean val- 
ues, can be studied in terms of the apparent standard devi- 
ation. The apparent standard deviation allows to estimate 
the mean deviations of characteristic dust quantities due 
to the turbulent fluid field. 

Equation H19() is therefore the mean, quadratic weighted 
deviation of the realisations i {i = O...Nx) of the tur- 
bulent, dust forming gas flow in space at an instant of 
time t. Figure depicts the space means (solid) {a)x{t), 
and the respective apparent standard deviations leading 
to {a)x{t) + <7^^_,it) (dotted) and (a),(i) - a^^_^{t) 
(dashed). Note that there is no straight forward functional 
dependence among all the lower (dashed) and all the up- 
per (dotted) curves, respectively. 

Figure El shows that the standard deviation is largest 
in the period of most active dust formation, i. e. between 
0.05 s and 3 s for the Mach number case depicted (compare 
paragr. "Dependence on Mach number"), independent on 
the initial reference temperature. For example, the mini- 
mum and the maximum deviations in density deviate by 
almost a factor 3 for the model with T^-d = 2100 K {y 
entree A in Table IX!2|l . The apparent standard deviation 
of the velocity field shows that Su{x, t) = ±0.2vrci and less 
which is subsonic (compare Figs. El 01 a-nd in agreement 
with Ludwig et al. (2002). 



The apparent standard deviations indicate that there 
are no very large deviations in the ID dust quantities if the 
dust complex has reached it steady state, in contrast to 
the hydrodynamic quantities. Turbulent fluctuations will 
cause the dust formation to set in somewhat earlier and 
to occur somewhat more vivid (larger J»). On the con- 
trary. Figs. El and El illustrates that no dust forms if the 
fluctuations result in T > T5 (no dashed line for Jn.). 



Dependence on temperature: The standard deviations of 
the hydrodynamic quantities are not considerably larger 
neither in a deeper nor in a shallower turbulence excited 
atmospheric layer (r.h.s. Fig. EJ. In contrast, the varia- 
tions in the dust quantities decrease with decreasing T^ef 
and increase with increasing Trof (the latter is not shown 
here). The nucleation rate decreases by orders of magni- 
tude with increasing temperature and the standard devia- 
tion is considerably larger. Consequently, the mean num- 
ber of dust particle is smaller and therefore the mean par- 
ticle size larger. In contrast, nucleation occurs earlier and 
more vivid with decreasing temperature. The dust forma- 
tion process is complete already during the very first time 
of the simulation as e. g. depicted on the r.h.s. in Fig. El 
therefore J, and Xhot are not resolved on the time interval 
plotted. Consequently, the number of dust particles is or- 
der of magnitudes higher and the mean particles size has 
decreased since the material has to be distributed over a 
larger grain surface area than if less particles form. All of it 
is in accordance with common expectations (see Sect . 13.411 . 
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Fig. 7. Time-means over r: 
tures Trei- dash-dot (green) 



lOmin for different tempera- 
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{reference results), dotted (red) - 1900K, dashed (blue) 
1500K (for details on reference values see Table ^ 



Dependence on Mach number; FigureElshows a simulation 
comparable to the r.h.s. of Fig. |S1 but now with M = 
0.02 instead of M = 0.1 (/ entrees A, B in Table rO|) . 
Consequently, the characteristic time scale is much longer, 
namely tref ~ 15 s instead of tref ~ 3 s and a much longer 
time interval needs to be depicted in Fig. El compared to 
Fig. [SI in order to observe the inset of the dust formation. 

Wc observe that the variation of the hydrodynamic 
quantities does not change remarkable compared to higher 
Mach number cases. It only appears on a much longer 
time scale. However, the superimposed waves need about 
3 times longer to initiate the first dust formation. The 
nucleation is somewhat less efficient resulting in a slightly 
lower number of dust particles which are, hence, slightly 
larger. Also the growth process is less efficient compared 
to the M=0.1-case 

Although the dust complex acts on its own, chemical 
time scales, it needs much longer time to reach a steady 
state situation if the initial Mach number is small (see also 
Fig. H in Sect. EH). 

3.4. The dust formation window 



the local temperature is already smaller than the nucle- 
ation threshold temperature. 

The stochastic regime contains those TD states for 
which the dust is possible if some realistic ignition mech- 
anism can cause T < Tg. This regime contains the critical 
range where a transition from T > Tg to T < Tg is pos- 
sible. The size of the stochastic regime depends on the 
turbulent energy. 

The third regime can be called impossible since no dust 
formation would be possible here. 

Temperature dependence: We have investigated the tran- 
sition deterministic ~ stochastic regime by studying the 
temperature dependence in our stochastic ID simulations. 
The time mean values (Fig. of the dimensionless hydro- 
dynamic variables are very much alike {{p)t, {p)t, {T)t, 
{u)t, for reference values see Table ^ but the dust quan- 
tities deviate considerably between these two extreme 
regimes^. 

Figure|3depicts four cases of which Tref = 1900, 1500K 
(dotted, dashed) fall into the deterministic regime in 
which the dust formation process is complete after a very 
short time on the whole test volume (/xi = fsi = 1) 
without any external excitation necessary. The Trcf = 
2500,2100K (dash-dot, sohd) faU into the stochastic 
regime where turbulence initiates the dust formation pro- 
cess by causing the very first nucleation event to occur 
(compare also Sect. P). In Tj-cf = 2100K-case the dust 
formation process is still completed after a very short 
time (t w 0.07s) in the whole test volume while for 
Tref = 2500K the first efficient nucleation event occurs 
only after about 65 s « Imin. The dust formation is not 
complete (/Ti,/si < 1) in this comparably hot case and 
even much more refined in time and space: The time mean 
of the mean particle size, ((a))t, varies by « 1 order of mag- 
nitude (4th panel, l.h.s. Fig.©. Therefore, T^f = 2500 K 
falls in the very end of the stochastic regime being close 
to impossible. 
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Stochastic fluctuation can drive a reactive gas flow into the 
dust formation window, i. e. the thermodynamic regime 
where the gas - solid (or liquid) phase transition is possible 
and most efficient (see e.g. Sedlmayr 1997). 

Depending on the thermodynamic (TD) situation, 
three regimes (compare Fig. Ill|l appear to be present in 
a turbulent atmosphere: The deterministic (or subcritic) 
contains those TD states where dust formation occurs 
without the need of an (e. g. hydrodynamic) ignition, i.e. 



The temperature sequence depicted in Fig. [3 displays 
the transition from the deterministic into the stochas- 
tic regime: the dust formation is most efficient at the 
smallest temperature considered (1500K) resulting in the 



The difference of {p)t in the case of T = 2500K is correct 
since here TnE/Tiai had to be considerably smaller in order to 
allow the system to enter the dust formation window (compare 
Table Q. 




Fig. 8. 2D simulations with Ti-cf = 2100 K, M = 1, TVk = 500( / entree A' Table rO|) for three instants of time 
during the period of active dust formation (top: 0.8s, middle: 1.7s, bottom: 8s; Ny^ x iVy ~ 128 x 128 = 500 m x 500 m, 
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most dust particles (3rd panel, r.h.s., Fig. O and there- 
fore in the smallest grain size. With increasing temper- 
ature, less particles are formed which can accumulate 
considerably more material and grow therefore to largest 
sizes. The hottest case considered seem not to fit into 
this picture because e.g. {{a))t,2500K > {{a))t,i500K but 
{{a))t,2500K ~ ((«))*, 1900K- Comparing the radiative equi- 
librium temperature of our test calculations (Table ^ in- 
dicates that Tre,250ok ~ 7re^i9ook- Apart from the fact 
that considerably more time is needed to form and grow 
the first, initiating seed particles the higher the tempera- 
ture is, the dust will drive the system locally towards this 
low radiative equilibrium temperature providing thereby 
TD conditions comparable to the T = 1900K-case. The 
resultant dust quantities need therefore be comparable if 
T = Tre, i. e. if the gas has reached the same isothermal 
state. 

3.5. 2D results 

So far, only ID results have been presented in this pa- 
per, which provide a good possibility to study the most 
important ongoing physical and chemical processes and 
their interactions. In ID, however, each wave crosses the 
whole test volume and will therewith influence the local 
thermodynamic conditions everywhere inside the volume. 
In 2D, the influence of waves leads to much more compli- 
cated patterns since e. g. a non-zero rotation of the fluid 
field can develop. The expected consequence is a much 
more heterogeneous distribution of dust than in any ID 
situation as it is illustrated in the following. 

A 2D model calculation with T^f = 2100 K, M = 1, 
A'^k^SOO ( y entree A' in Table IA.2p was performed on 
a spatial grid of x Ny = 128 x 128 cells correspond- 
ing to a box of 500m x 500m. The smallest eddies have 
a size of A^^^ = 5 m, the largest are of the size of the 
test volume. The gravity acts in negative y-direction, i. e. 
g = {0,-3,0}. The initially homogeneous and dust free 
fluid is constantly disturbed by superimposed waves enter- 
ing from the left, the right, and the bottom side. Thereby, 
a gas element is modelled which is continuously disturbed 
by waves originating from the surrounding convectively 
instable atmospheric fluid. The top side is kept open sim- 
ulating the open upper boundary of a test volume in the 
substellar atmosphere. 

Figure IHl shows three instants of time during the phase 
of vivid dust formation and demonstrates the appear- 
ance of large and small dusty scale structures evolving 

^ Low Mach number simulations of driven turbulence in 2D 
are not yet possible with the present code. Botta et al. (2003) 
have shown that unbalanced truncation errors can lead to con- 
siderable instabilities in the complete, time-dependent equa- 
tion of motion in a quasi-static situation and suggest a bal- 
anced discretisation scheme. We will tackle this problem in a 
forthcoming paper and use our present 2D results for M=l 
only to illustrate the stronger influence of the hydrodynamic 
processes on the evolving dust structures to be expected in 
multi-dimensional simulations compared to ID. 



with time. Both, the number of dust particles Ud (l.h.s.) 
and the mean particle size (a) (r.h.s.), are plotted on a 
logarithmic scale with = 1 . . . 10^ cm~^ and (a) = 
10~^-^ . . . 10~^'^. The very inhomogeneous appearance of 
the dust complex is a result of nucleation fronts and nucle- 
ation events comparable to our ID results. The nucleation 
is now triggered by the interaction of eddies coming from 
different directions. Large amounts of dust are formed and 
appear to be present in lane-like structures (large logna; 
dark areas). The lanes are shaped by the constantly in- 
ward travelling waves. Our simulations show that some 
of the small-scale structure merge thereby supporting the 
formation of lanes and later on even larger structures. The 
formation of such large structures is not caused by the es- 
tablishment of a pressure gradient to counterbalance the 
gravity. Since the whole test volume is only of the size of 
Hp/ 20 the resulting pressure gradient is negligibly small. 
Hence, the large-scale structures result from the interac- 
tion of dust formation and turbulence. 

Furthermore, dust is also present in curl-like structure 
which indicates the formation of vortices. As the time 
proceeds in our 2D simulation, vortices develop orthog- 
onal to the velocity field which show a higher vorticity 
(V X v{x,t)) than the majority of the background fluid 
field. For illustration, the maximum and the minimum vor- 
ticity between w — 20s~^ and ~ 20 s~^ has been superim- 
posed as a contour plot (grey/black) on top of the false 
colour plot of the number of dust particles for t = 0.8 s 
in Fig. |S1 (l.h.s., top). This shows that the vortices with 
high vorticity preferentially occur in dust free regions or 
regions with only little amounts of dust present. The mo- 
tion of the vortices can transport the dust particles into 
region with still condensible material available, and seem 
thereby to cause larger and larger dusty areas to form. 

The 2D hydrodynamic behaviour is comparable to our 
ID results in the sense that the stochastically created 
waves enter the test volume, interact, run through the 
test volume, and eventually leave it at the top side af- 
ter they have initiated the dust formation process. Large 
changes in the dust quantities during short time inter- 
vals occur locally during the time period of the first nu- 
cleation and the re-initiation of nucleation by radiative 
cooling. A spatially inhomogeneous dust distribution re- 
sults which is only slightly shifted back and forth by the 
inward moving superimposed waves. In contrast, consid- 
erable variations in p, p, and in the velocity components 
(u, v) occur in time and space. Furthermore, a new quali- 
tative behaviour compared to the ID simulation appears. 
Hydrodynamic advection gathers the dust in larger and 
larger structures while the formation of the dust has been 
initiated in the small-scale regime of the simulation. The 
reference time of our present 2D simulation is compara- 
bly high (see Table IX!2|l . We have, however, done ID test 
calculations (Figs. [Sj El EJ which indicate that in the low 
Mach number case the global dust formation time scale 
(i. e. the time when the dust formation reaches its steady 
state) will merely increase (for discussion see Sect. I4.2|l . 
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Fig. 9. Space mean and standard deviations (solid) of the total grey absorption coefRcient k = Kgas + Kdust and the 
total grey optical depth t for the simulations depicted in Fig. The dust (dashed) and the gas (dotted) contributions 
are shown for simulations with Tj-ef = 2100 K with M = 0.1 (l.h.s.; entree A Table IA.2(l and with M = 0.02 (r.h.s.; 
entree B Table [0)) . 



4. Discussion 

4.1. Variability 

Observational evidence has been provided (Bailer- Jones 
& Mund 1999, 2001a,b; Martin et al. 2001; Gehno et 
al. 2002; Clarke et al. 2002) that Brown L-Dwarfs are 
non-periodically photometric variable at a low level of 
1 — 2% (Clarke 2002) and sometimes even larger. Nakajima 
et al (2000) and Kirkpatrick et al. (2001) reported on 
spectroscopic variability. Appealing explanations are the 
appearance of magnetic spots or the formation of dust 
clouds. Mohanty et al. (2002) and Gelino et al. (2002) 
have argued that ultra cool dwarfs are unlikely to support 
magnetic spots. There is, however, evidence for magnetic 
activity in L dwarfs by a rapid declining, strong Ha emis- 
sion. In contrast, three objects are observed with a persis- 
tent, strong Ha emission (Licbert et al. 2003). A straight 
forward, consistent explanation is not at hand yet and will 
probably be theoretically very demanding. We, therefore, 
follow in this paper the hypothesis of the formation of dust 
clouds in a convcctively influenced turbulent environment 
as explanation of non-periodic variability 

Ludwig et al. (2002) argue based on their hydro- 
dynamic 3D simulation that M-dwarfs (and even more 
Brown Dwarfs) show only very little temporal and hor- 
izontal fluctuations in their atmospheres. Dust strongly 
interacts with the thermo- and hydrodynamics due to ra- 
diative transfer effects, gas phase depletion and on macro- 
scopic scales due to drift. It seems therefore likely to 
expect a support of initially small inhomogcneouties by 
these processes. Woitke (2001) has carried out 2D radia- 



tive transfer calculations for a inhomogencous density dis- 
tribution which support this idea. Since the radiation is 
blocked by condensing dust clouds of sufficient optical 
depth, the radiation is forced to escape mainly through 
the remaining wholes, thereby enhancing and preventing 
the dust formation, respectively. 

One may speculate that the consideration of dust for- 
mation in 3D simulation may even cause these models to 
deviate considerably from the simple MLT models in the 
case of brown dwarfs due to the time dependence of the 
dust formation process and the corresponding feedback on 
the space and time evolution. 

Optical depth: Figure |51 (left panels) depicts the space 
mean total opacity (solid lines) and the mean opacity of 
the dust (dashed lines) and the gas (dotted lines) for the 
ID test calculation investigated in Fig. |7| Upper curves 
indicate {K)x{t) + cr^^_]^(t), the lower depict {K)x{t) — 
''^Nt-i(^)- "^^"^ Rosseland mean dust opacities for astro- 
nomical silicates (Kdust = 0.75 pL^l.lAT^-^^, Paper I; 
l.h.s.) and a typical Rosseland gas mean of Hg^s = 

0. 1g/cm'^ have been adopted. The Rosseland gas mean 
opacity was chosen typical for a hot, inner layer of a brown 
dwarf atmosphere. Figure shows that the dust and the 
gas opacities differ by about 1.5 orders of magnitude. The 
dust opacities varies by about 0.5 magnitudes, the gas 
opacity by only about 0.2 order of magnitudes, this is in- 
dependent of the characteristic time scale of the system, 

1. e. the large-scale Mach number of the initial configura- 
tion (compare l.h.s. and r.h.s. in Fig.l^j). 
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Figure El further depicts in the right panels the time 
evolution of the space mean of the optical depth which a 
mass element of the size of our test volume (500 m) would 
have. The optical depth increases by a factor of 10 when 
the dust forms in accordance with the opacity increase. 
There is, however, a delay between the inset of dust for- 
mation and the time when the maximum optical depth 
is reached because the dust formation sets in earlier than 
a cloud becomes optically thick. This delay is sensitive to 
the characteristic time scale of the system as a comparison 
of the l.h.s. and r.h.s. in Fig. [^demonstrate. 

The gas/dust mixture will only be optically thick in 
the case of non-transparent dust particles (e. g. astronom- 
ical silicates) which may depend on the wavelength con- 
sidered. Glassy grains may efficiently absorb in the far IR 
(A <; 10/im) and one might, depending on the observed 
wavelength, detect typical dust features and possibly even 
the rapid formation events which, however, may need to 
evolve on macroscopic scales in order to be observable. 

Turbulent fluctuation seem not to produce consider- 
able variations in view of the space mean of the optical 
depth in the long term (here t > As) if astronomical sili- 
cates are assumed as typical dust opacity carriers. 

4.2. Towards the observable regime 

The time-dependent simulations performed so far might 
suggest an up-scaling of the achieved results in order to es- 
timate e.g. possible cloud sizes and variability time scales 
on an observational level. This idea is not as straight for- 
ward as it might seem because of the differences between 
chemical and hydrodynamical time scales. Dust formation 
is a local process and its (local) time scale is the same in- 
dependent of the hydrodynamical regime but the hydro- 
dynamic time scale changes depending on e. g. the char- 
acteristic length of the regime considered. A simple and 
correct up-scaling according to the principle of similarity 
would require the characteristic numbers to remain the 
same which is not the case (see also Helling 2003, Lingnau 
2004), for instance due to the changing characteristic hy- 
drodynamic time scales. One might further imagine to up- 
scale the results by a periodic continuation, i.e. each small- 
scale simulation being one pixel of the large-scale picture. 
In any case, the feedback of the large scale (e.g. convec- 
tive) motions on the small scales is neglected and visa 
versa. Important effects like chemical mixing and kinetic 
energy input into the turbulent fluid field are than miss- 
ing. Consequently, the non-linear coupling between dust 
formation and hydrodynamics makes it impossible to sim- 
ply scale up detailed small-scale results into an regime 
accessible by observations. The complete large-scale sim- 
ulation needs to be awaited. 

The aim of the investigations of the small-scale regimes 
performed so far was indeed to provide information and 
understanding for building a large-scale model of a brown 
dwarf atmosphere. From these results, we suggest the fol- 
lowing necessary criteria for a sub-grid model (also called 
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Fig. 10. The deterministic and the stochastic dust for- 
mation regime in comparison to typical substellar model 
atmospheres (solar metalicity, \ogg = 5, AH (black) = 
AUard et al. 2001, T (grey/cyan) = Tsuji 2002). 
(open circles (black) - dust formation without ignition; 
asterisk (blue) - dust growth on initially present seed par- 
ticles (Icm^"^); squares (red) - dust formation under tur- 
bulent conditions; triangles (black) - (Trcf , p) pairs inves- 
tigated; hooks - border of convcctive zone) 



closure term or closure approximation) of a turbulent, dust 
forming system: 

a) A sub-grid model must describe the transition deter- 
ministic —^ stochastic dust formation, depending on 
the turbulent energy as e. g. measured by the large- 
scale variations of the local velocity field. 

b) The dust formation process (nucleation -I- growth) is 
restricted to a short time interval (of the order of a few 
seconds) , which is usually much smaller than the large- 
scale hydrodynamic time scale. This involves that: 

• The nucleation occurs locally and event-like in very 
narrow time slots. 

• The growth process continues as long as conden- 
sible material is available and thermal stability of 
the dust is assured. 

• The condensation process finally freezes in and the 
inhomogeneous dust properties are preserved. 

c) The dust formation process should be accompanied by 
a fast transition from an approximately adiabatic to 
an approximately isothermal behaviour of the dust / gas 
mixture. This transition can be expected to affect the 
convective stability in substellar atmospheres. 
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4.3. Comparison with classical model atmospheres 

Figure [Tni revisits the idea of the dust fornration window 
but now with view on classical brown dwarf atmosphere 
calculations. 

The following calculations are performed for various 
(T, p) pairs in order to demonstrate in which atmospheric 
regions dust formation will simply take place (determin- 
istic regime) and in which regions dust formation needs 
to be initiated (stochastic regime; compare Sect. l3.4f) . For 
eye guidance, the (T, p) pairs are indicated by small black 
triangles along with presently used substellar model atmo- 
spheres in the literature (AUard et al. 2001 - black lines, 
Tsuji 2002 - grey/cyan lines). One may notice that the 
warm models of the same stellar parameter can differ by 
about 500 K for a given density in the inner atmospheric 
region which is convectively unstable in both cases. Small 
hooks (red) indicate the upper Schwarzschild boundary 
(SB) for convection if available, i. e. wmlt = for T < T^b- 
The test is simple: Two calculations are carried our for 
each (T, p) pair where in the first test a small number 
of seed particles is prescribed {pLg = Icm"'^; open black 
circles), and in the second test no seed particles are pre- 
scribed (blue asterisks) . The results should be almost iden- 
tical inside the deterministic regime but no dust should 
form in the stochastic regime. 

One observes that there is no difference in grain size in 
the cool outer atmospheric region whether seed particles 
are initially present or not. The initially present little num- 
ber of seed particles is overran by a very efficient nucle- 
ation which is characteristic for the deterministic regime. 
Note that the dust grains will immediately start to move 
inwards on macroscopic scales of atmospheric extension 
due to the large gravity of brown dwarfs. However, de- 
pending on the TD conditions, the nucleation efficiency 
varies in the atmosphere. The more efficient nucleation 
takes place, the more seed particles are formed resulting 
in smaller mean sizes because the gaseous material has 
to be distributed over more particles than in case of less 
efficient nucleation. Small mean particle sizes inside the 
deterministic regime are therefore a sign for efficient dust 
nucleation. 

For T > 1800 K, the stochastic regime is entered. Here, 
dust is only present if by some mean a certain number of 
seeds is present (blue asterisks) or a ignition mechanism 
(e. g. turbulence, radiative cooling by gas) provides the 
appropriate TD conditions (red squares). Note that no 
open circles appears at this site of Fig. ^| since the TD 
conditions are now inappropriate for a gas - solid phase 
transition except an ignition mechanism is present. 

5. Conclusions 
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Fig. 11. Regimes of turbulent dust formation. T5: nucle- 
ation threshold temperature (supersaturation 5 ^ 1 re- 
quired), Tsub: sublimation temperature. 



lent motions and just reaches sufRciently low temperatures 
for condensation, but other applications are also conceiv- 
able. The dust formation in a turbulent gas is found to 
be strongly influenced by the existence of a nucleation 
threshold temperature Tg. The local temperature T must 
at least temporarily decrease below this threshold in order 
to provide the necessary supersaturation for nucleation. 

Depending on the relation between the local mean tem- 
perature T and Ts, three different regimes can be distin- 
guished (see l.h.s. of Fig.[TT|l: (i) the deterministic legiine 
(T < Ts) where dust forms anyway, (ii) the stochastic 
regime (T > Ts) where T < Ts can only be achieved 
locally and temporarily by turbulent temperature fluctu- 
ations, and (iii) a regime where dust formation is impos- 
sible. The size of the stochastic regime depends on the 
available turbulent energy. This picture of turbulent dust 
formation is quite different from the usually applied ther- 
modynamical picture (r.h.s. of Fig. Ill|l where dust is sim- 
ply assumed to be present whenever T < Tguh, where Tgub 
is the sublimation temperature of a considered dust ma- 
terial. 

The investigations have been performed in the meso- 
scopic scale regime, where the test volume is excited by a 
spectrum of waves within a limited fc-interval at given en- 
ergy distribution, which are generated at the boundaries 
(pseudo-spectral method for driven turbulence). Two ba- 
sic processes are found to be capable to initiate the dust 
formation even in a host-hostile environment: 



We have studied the onset of the dust formation process 
in a turbulent fluid field, typical for the dense and ini- 
tially dust-hostile regions in substellar atmospheres. The 
main scenario is a convectively ascending fluid element 
in a brown dwarf atmosphere, which is excited by turbu- 



1) expansion waves with T— AT < Ts {nucleation fronts) . 

2) interactions of two or more expansion waves which can- 
not produce sufficiently low temperatures for them- 
selves, but the superposition of such waves can T — 
ATi — AT2 — . . . < Ts {nucleation events). 
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After initiation, the dust condensation process is com- 
pleted by a phase of active particle growth until the con- 
densible elements are consumed, thereby preserving the 
dust particle number density for long times. However, ra- 
diative cooling (as follow-up effect) is found to have an 
important influence on the subsequent dust formation, if 
the dust opacity reaches a certain critical value. This cool- 
ing leads to a decrease of T{t) which may re-initiate the 
nuclcation. This results in a runaway process (unstable 
feedback loop) until radiative and phase equilibrium is 
achieved. Depending on the difference between the ini- 
tial mean temperature T(t = 0) and the radiative equilib- 
rium temperature Tre, a considerable local temperature 
decrease and density increase occurs. Since the turbulent 
initiation of the dust formation process is time-dependent 
and spatially inhomogcncous, considerable spatial varia- 
tions of all physical quantities (hydro-, thermodynamics, 
dust) occur during the short time interval of active dust 
formation (typically a few seconds after initiation), which 
actually creates new turbulence. 

Thus, small turbulent perturbations have large effects 
in dust forming systems. A convectively ascending, ini- 
tially dust free gas element, which is slightly warmer than 
its surroundings, can be excited to form dust by waves 
running through it, even at otherwise dust-hostile tem- 
peratures. The newly created dust particles may cause the 
substellar atmosphere to become almost instantaneously 
optically thick. Our 2D simulations show that the dust ap- 
pears in lane- like and curled structures. Small-scale dust 
structures merge and form larger structures. Vortices ap- 
pear to be present preferentially in regions without or with 
only little dust. Non of these structures would occur with- 
out turbulent excitation. 
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Appendix A: Analysis of characteristic numbers & 
characteristic values of test calculations 
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Table A.l. Characteristic numbers and reference values used to analyse the chemistry and physics in the mesoscopic 
scale regime. Hydro- and thermodynamic reference values have been guided by static brown dwarf model atmosphere 
calculations (Tsuji 2002). The reference values for the dust complex have been adopted from the experiences gained 
in Paper I. Considered is the Ti02 seed formation and growth by SiO and Ti02. (7 = 7/5, aj^^'' = 1.95 ■ 10~^cm, 
mu mass of the hydrogen atom, a Stefan-Boltzmann constant, k Boltzmann constant, /^kin kinematic viscosity [g 
cm-^s-i], Ni = 1000) 



Name 


Characteristic 
Number 


Value 

inside outside 


Reynolds number 

Mach number 
Froude number 


Rc - 

Mkin 

K 00 in — ^ pT^ (*) 
— 0.08 • iU V-t 

M= 
Fr = 


101 aIO 101 n"^ 
i.z iU ... i.z iU 

0.1 

0.11 0.05 


Radiation number 




n 907 SRin~3 (<>) 

U.ZU ( ... o.D lU 

2.07 ... 0.86 


Damkohler number of nucleation 
Damkohler number of growth 

Sedlmayr number {j £ N) 


n„gr _ iiofXrof 


... 3.2410^ 

9.78 lO"* ... 97.8 

i = : 1 
i = 1 ; 0.195 
3 = 2: 0.0381 
j = 3 : 7.44 10"^ 


FlpTTipnt r^nriQi] mntinn vniTTibpT" 


_ Prof -i^Crof 

"<H>.rct<:rct 


7.27 10""* 


Name 


Reference Value 


Value 

inside outside 


temperature 
density 
thermal pressure 

velocity of sound 
velocity 
length 
hydrodyn. time 
gravitational acceleration 


r,ef [K] 

Pref [g/cm^] 

^^cf = [dyn/cm^] 
cs = -y7^ [cm/s] 

Href [cm/s] 
Zref [cm] 
tref = ^ [S] 
flref [Cm/S^] 


2200 . . . 1000 
10~=* ... lO"'' 
7.7810^ ... 3.5410* 

3.310^ ... 2.2310^ 

^cs/10 
10^ 

3.03 . . . 4.49 

lO'"^ 


total absorption coefRcient 


«rof [1/cm] 


10"^ ... 3-10~'^ 
0.1 ... 3-10"^ 


nucleation rate 
O''' dust moment {— rid/ p) 
heterog. growth velocity 
mean particle radius 


J*,rof/n<H>,ref [1/s] 
i0,rcf [1/g] 

Xrcf [cm/s] 
{a)rcf [cm] 


... 2.5 -lO"" 
1.35 10^* 

<0 ... 3.5110"^ 

10"'' 


element abundance 
total hydrogen density 
SiO density 
Ti02 density 


Ercf [-] 

n<H>.rof = [i/cm^] 

w<SiO>,rcf [1/cm^] 

n<Ti02>,rcf [1/cm^] 


10"" 

4.22 10^° . . . 4.22 10^^ 
1.50 10^'' ... 1.5010" 
4.1210" ... 4.1010" 



For a H2/He rich gas (Eq. (10) in Woitke & Helling 2003a). gas only dust and gas Gail & Sedlmayr (1999) 
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Table A. 2. Characteristic numbers and reference values for the test calculations carried out in the previous sections. 
For definition of characteristic numbers see Table 





A 


B 


c 


A' 






Hydro- & Thermodynamic 






Trof 


[K] 


2100 


2100 


1900 


2100 


Prcf 


[g cm~^] 


3.1610-" 


3.16 10"* 


8.1310-" 


3.16 10"* 




[cm] 


lO'^ 


lO'^ 


10^ 


10^ 


"rcf 


[cm s""""] 


3.25 lO"* 


6.50 10^ 


3.09 10" 


3.25 lO'' 


9rcf 


[cm s^^] 


10^ 


10^ 


10^ 


10^ 


The 


[K] 


1980 (= o.gr.ef) 

210"^ 


1980 (= 0.9rrcf) 


1710 (= 0.9T,rf) 


1980 (= 0.9rrcf) 


^rcf 


[1/cm] 


210"^ 


210^^ 


210"^ 






Dust & Chemistry 






r 

* ,rcf 


-3 

[s cm ] 


10^2 


10^^ 


10^2 


10^2 


Iref 




10"^ 


10"® 


1 n"*^ 


1 


Aref 


[cm s ^] 


10"^ 


10"^ 


10"^ 


10"^ 


T „ r 

J^O,rci 


IS J 


10^" 


10"^ 


10"' 


10^° 






Derived reference values 






Prci 


[dyn cni~'^] 


2.38 10^ 


2.38 10^ 


5.54 lO'' 


2.38 10^ 


f ^ 

'TCI 




3.08 


15.4 


3.24 


0.3 


ij.rcf 


r i — li 

[cm-ig ] 


( 3"=' )■" -Lo.rof 




( 3^) ' io.rcf 


/ 47r f \773 r 

( Lo.rcf 






Hydro- & Thermodynamic 






M 




1/10 


1/50 


1/10 


1 


Fr 




0.105 


4.22 10"^ 


0.095 


10.5 


Rd 




1.139 10"^ 


5.70 10"^ 


3.45 10"^ 


1.1410"^ 






Dust & Chemistry 






DaT" 




9.73 10'^ 


4.8710'' 


3.9810^ 


9.73 10" 


Daf 




1.91 lO"* 


9.55 lO'' 


2.01 10" 


1.91 10^ 


Seo 




1.0 


1.0 


1.0 
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Sei 




0.195 


0.195 


0.195 


0.195 


Se2 




3.810"^ 


3.810"^ 
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7.41 10"^ 


7.41 10"^ 


7.41 10"^ 


7.41 10"^ 


El 




1.41 10"'' 


7.0710^ 


3.82 lO"'' 


1.41 10"^ 






Turbulence 






e (Eq.EJ 


[cm^ s"^] 


2.40 10** 


1.92 lO'' 


2.06 10** 


2.41 10" 


Srof (Eq.Ei 


[erg/K] 


1.42 10^ 


1.3710^ 


1.2810^ 


1.410^ 


^min — 27r//max 


[1/cm] 
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fcmax = 2TT/{3dx) 


[1/cm] 
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[cm] 


^rcf 


^rcf 


^rcf 


^rcf 
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[cm] 
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dx for iV^ = 128 (2D) 


[cm] 








3.9410^ 
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